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Abstract. Despite originating in science fiction, warp drives have a concrete 
description in general relativity, with Alcubierre first proposing a spacetime metric 
that supported faster-than-light travel. Whilst there are numerous practical barriers 
to their implementation in real life, including a requirement for negative energy, 
computationally, one can simulate their evolution in time given an equation of state 
describing the matter. In this work, we study the signatures arising from a warp drive 
‘containment failure’, assuming a stiff equation of state for the fluid. We compute the 
emitted gravitational-wave signal and track the energy fluxes of the fluid. Apart from 
its rather speculative application to the search for extraterrestrial life in gravitational- 
wave detector data, this work is interesting as a study of the dynamical evolution and 
stability of spacetimes that violate the null energy condition. Our work highlights the 
importance of exploring strange new spacetimes, to (boldly) simulate what no one has 
seen before. 
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1. Introduction 


Detections of gravitational waves by the LIGO Scientific-Virgo-Kagra Collaborations, 
e.g., [1-3] are paving the way to viewing extreme gravitational events across the entire 
visible universe. It is natural to ask what new signals originating from strongly distorted 
regions of spacetime could be seen in the future beyond the compact binary mergers 
already detected. To search for such signals and to correctly identify them in the 
measured data, it is important to understand their phenomenology and properties. 
Numerical relativity (NR) is one of the main tools for computing gravitational-wave 
signals for spacetimes in the highly dynamical strongly gravity regime in the absence of 
symmetries [4-7]. NR simulations are time domain evolutions of a given initial state that 
allow the resulting gravitational-wave signals to be extracted in the asymptotic radiation 
zone, from which we can connect to our own observations from a given viewing angle 
and distance. 

Unlike black holes, which have a natural formation channel [8, 9] and considerable 
observational support [10-12], more exotic spacetimes such as wormholes or warp drives 
[13, 14] are considered science fiction because their formation and existence generate 
many potential paradoxes, and require matter that violates certain energy conditions. 
Nevertheless, warp drives can be described classically within the theory of general 
relativity by a metric first proposed by Miguel Alcubierre [15] and further developed in 
subsequent works, e.g., [16-19]. 

The principle idea behind a warp drive is that instead of exceeding the speed of 
light directly in a local reference frame, which would violate Lorentz invariance, a “warp 
bubble” could traverse distances faster than the speed of light (as measured by some 
distant observer) by contracting spacetime in front of it and expanding spacetime behind 
it [20]. The problems associated with such spacetimes are well documented [13, 21- 
25]. One particular issue is the requirement for matter that violates the Null Energy 
Condition (NEC) [26, 27] {. The requirement that warp drives violate the NEC may be 
considered a practical rather than fundamental barrier to their construction since NEC 
violation can be achieved by quantum effects and effective descriptions of modifications 
to gravity [32, 33], albeit subject to quantum inequality bounds [34-36] and other 
semiclassical considerations [37] that seem likely to prove problematic. Other issues 
with the warp drive metric include the potential for closed time-like curves [21, 23] and, 
from a more practical perspective, the difficulties for those in the ship in controlling and 
deactivating the bubble [38]. 

From the perspective of simulating the warp drive dynamically, the key challenge 
is stability. One can straightforwardly generate the initial metric that describes the 
Alcubierre metric and infer the instantaneous matter configuration that must support it 
using the Einstein Equation, but a time-domain evolution of the coupled matter-gravity 
system necessitates that we specify the equation of state of the matter, which in turn 


t Recent works that have suggested positive energy warp drive solutions [28-30] have problems that 
are discussed in [31]. 
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determines how it evolves. In a perfect fluid, this involves specifying the relation between 
the fluid pressure and density, but in a more general case, the full relation between the 
stress-energy components must be given. There is (to our knowledge) no known equation 
of state that would maintain the warp drive metric in a stable configuration over time 
- therefore, whilst one can require that initially, the warp bubble is constant, it will 
quickly evolve away from that state and, in most cases, the warp fluid and spacetime 
deformations will disperse or collapse into a central point. 

This instability, whilst undesirable for the warp ship’s occupants, gives rise to 
the possibility of generating gravitational waves. The Alcubierre metric for a constant 
velocity warp ship would not generate gravitational waves. In addition, far away from 
the warp bubble, the spacetime is exactly Minkowskian — that is, it is actually flat, not 
just asymptotically flat. The warp bubble spacetime is therefore undetectable by its 
gravitational effects as it has no gravitational wave content and no measurable ADM 
mass [39] §. However, the collapse of the warp drive, or its acceleration or deceleration, 
should generate gravitational waves, and in this work, we study the first case by allowing 
the warp bubble to collapse. Physically, this could be related to a breakdown in the 
containment field that the post-warp civilisation (presumably) uses to support the warp 
bubble against collapse. 

The prospect of a zero ADM mass spacetime radiating gravitational-wave energy 
immediately leads to an interesting physical question: Does the spacetime actually 
radiate energy, and if so, does the final spacetime have a negative ADM mass? In this 
work we consider the fluxes of energy both from the matter and the gravitational waves 
out of a spatial volume after the warp drive collapse. This study has connections to 
early work by Bondi [42], in which positive and negative matter restricted to linear 
motion was shown to attract and repel in a way that would accelerate them together. 

For our study, we restrict ourselves to a limited and simplified scenario and sub- 
light speeds as a proof of principle||. The initial bubble is described by the original 
Alcubierre metric with a fixed wall thickness. We develop a formalism to describe the 
warp fluid and its evolution and vary its initial velocity at the point of collapse (which 
is related to the amplitude of the warp bubble). We neglect the existence of the warp 
ship itself [27], which is assumed to be negligible compared to the bubble. 

This article is structured as follows. In Sec. 2, we set out the formalism used 
to perform our simulations of warp bubble collapse, particularly the treatment of the 
warp fluid element and its equation of state. In Sec. 3, we present the resulting 
gravitational-wave signatures and quantify the radiation of energy from the spacetime 
region. In Sec. 4, we discuss the implications of our findings and propose directions for 
further study. Some technical details and validation of our simulations are given in the 


§ It circumvents theorems about all asymptotically Minkowskian spacetimes being Minkowski [40, 41] 
by the presence of the negative energy contribution. 

|| Our set-up appears to work for speeds of order v < 1 but requires much higher resolutions for larger 
speeds. It would be better to use a code adapted to the cylindrical symmetry of the problem to study 
such cases. 
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Appendix. 
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Figure 1. The evolution in time of the trace of K;ij, K = y" Kij. This quantity 
initially provides a schematic representation of the warp drive metric, with positive 
values for the contracting part of the space ahead of the ship and negative values 
for the expanding part behind. We observe the destabilisation of these regions, 
which initially collapse inwards and then dissipate to infinity. Movie at https: 
//youtu.be/6aqkxooXAc4?si=NhremdNv5_F3sTKL. 


2. Formalism for the warp drive simulation 


We follow the standard numerical-relativity approach to solving the Einstein Equation 
1 
Cw = Ry = 5 Iw = STT uv (1) 


as an initial value problem, using the 3+1-decomposed form and setting G = c = 1. This 
means that we define spacelike hypersurfaces that foliate the four-dimensional manifold, 
specify initial data on a particular hypersurface, and evolve it forward in our chosen time 
coordinate using the method of lines. (See e.g. [5, 43, 44] for detailed discussions.) The 
3 + 1-decomposed metric in coordinates adapted to the slicing takes the form 


ds? = (—a? + B,6")dt? + 28,dtda' + yidx'de’, (2) 


where a is the lapse function, 3" is the shift vector, and ¥;; is the spatial metric connected 
to the spatial hypersurface. The description is completed by the extrinsic curvature K;;, 
which is roughly related to the time derivative of y;;, and more rigorously defined as its 
Lie derivative, such that 


1 
Ki = -zza + DeBi- (3) 


It will also be useful to note that n, = (—a,0,0,0) is the normal to the spatial 
hypersurface. 


What no one has seen before: gravitational waveforms from warp drive collapse 5 


2.1. Initial data for the metric 


Although, as discussed above, many alternatives with improved physical properties have 
been discussed in the literature, we will consider in this work the original Alcubierre 
metric [15] for which 


a =], (4) 
Jij = ij, (5) 
B = (0, 0, =, ] (Tx) ) (6) 
with d 
wat, ma Verte tea), (7) 
and 
Fe eS aie) (8) 


The free parameters R and o determine the shape of the warp bubble, specifically its 
radius and the wall thickness ~ 1/ø. 

We will use the value of R to define the system of units in which our results are 
expressed - our simulations have a scaling freedom that means that they apply to any 
similar warp drive with size R. Specifying the physical value of R determines the other 
physical values of our measurements, including the gravitational-wave fluxes and the 
energy radiated. For reference, for a warp bubble of size R = 1 km, the relevant 
quantities in physical units are t = 1[R] is 3.33 us and E or M = 1[R] is 0.677 Mo $. 

The extrinsic curvature carries the information about the expansion and contraction 
of the spacetime, with values 


1 
Ki =5 (dib; + 0;8i) (9) 


that can be computed analytically from the expression for f(r,). The trace of K;ij, 
K = 7"K;; is initially K = 0,67. This quantity is usually plotted to provide a schematic 
representation of the warp drive, see Fig. 1. 


2.2. Initial data for matter 


The initial data for the matter degrees of freedom is computed from the energy- 
momentum-tensor that supports the initial spacetime curvature deformation. We 
describe the state of the matter at any time in terms of projections of the stress-energy 


tensor 
Spv = Eo oii (10) 
So. == Pon Ty, (11) 
P = nën” Tuv, (12) 


€| For comparison, the Enterprise-E of the Star Trek universe is 685.7 meters in length [45]. 
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Figure 2. The evolution in time of the real part of the Weyl scalar Y4 for the case 
of v = 0.1. This quantity provides a measure of the spacetime curvature and, in the 
far region, can be identified with the gravitational wave content of the spacetime. 
In the last two panels, we see a burst of gravitational-wave radiation leaving the 
collapsed remnant of the warp bubble. Movie at https://youtu.be/h-BieIMO9Xw? 
si=vVFQf0k9nZiQGQkz. 


with the projection operator P” = 64 + nn, where n” is the normal vector to the 
spatial hypersurface defined above. In coordinates adapted to the slicing, we can work 
with the ten components 


T= Ta (13) 
1 BI 
Sy = — a it + a (14) 
1 Bt oe 
Pp = gatt = ait + a Li (15) 


Knowing that the energy-momentum-tensor has to satisfy the Einstein Equation, we 
can directly calculate the quantities on the right-hand side from Eq. (1) as 
1 R 
Lig = = (ru — Fou) , (16) 
where the Ricci tensor can be straightforwardly calculated from the Alcubierre metric 
in terms of derivatives of the shift. 


2.8. Metric evolution 


For the metric evolution, we use the standard CCZ4 formulation of [46] and the moving 
puncture gauge [47-50] within the numerical-relativity code GRCHOMBO [51-53], which 
uses the method of lines, with an RK4 time integration and 4th order finite difference 
stencils for calculating spatial gradients. 

We note that for the gauge evolution we use the shock-avoiding Bona-Masso type 
gauge slicing originally proposed by Alcubierre [54, 55], and recently advocated for by 
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Figure 3. The evolution in time of the matter energy density p for the case of 
v = 0.1. We see that a ring of positive energy forms within the initially negative 
energy density bubble and propels it outwards, which is reminiscent of [42]. There 
are several waves of alternating positive and negative energies. In the final panel, 
we see that the remnant appears stable, and compared to Fig. 2, we see that the 
matter waves propagate at roughly the same speed as the gravitational waves. Movie 
at https: //youtu.be/RGRB5Y0c17w?si=e6z1n3_ijxRDNd44. 


Baumgarte and Hilditch [56] in the context of critical collapse simulations. We find that 
the shock-avoiding slicing gives much better stability and, more importantly, allows us 
to extract coherent gravitational waveforms in the wave zone, whereas in the standard 
gauge, we obtained large unphysical gauge effects that hindered us in extracting the 
physical gravitational-wave signal. Because this gauge has a tendency to push the 
extrinsic curvature to zero, one should not over-interpret the evolution shown in Fig. 1. 
It quickly decays to zero in the simulation before the matter and gravitational-wave 
fluxes have begun to leave the volume. 

Better quantities to track the evolution of the metric curvature are the Wey] scalars, 
in particular Y4. We study the tensorial gravitational-wave modes emitted by the 
warp drive by extracting the Newman-Penrose scalar Y, with the tetrads proposed by 
[57], projected into spin-weight —2 spherical harmonics, Wm = $52 Valr=rex [2Y] d9, 
where dQ = sin @ dé dy is the area element on the S? unit sphere. The evolution of the 
real part of W, is illustrated in Fig. 2. 


2.4. Matter evolution 


The requirement that V,7T“” = 0 gives us four evolution equations for the components 
of the matter p and S; as follows: 


; 1 z 
op = B Dip = Tan (a° Si) + a(S K” + pK), (17) 


OS; = BI DS; + SDBI — 7" Dla Sij) + aK S; — pdja. (18) 
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The evolution equation for S;; is, in principle, free - once we specify the equation of 
state of the matter, it would be fully determined. We base our choice on the observation 
that the initial warp drive spacetime has an average pressure of [31] 


a stig. (19) 


p= yË S; /3 = p — 


which we can reexpress using the ADM formalism as 


o p 1 /D’a z 
p-f- a Hyr” + K*) ine 


We then decompose the stress density tensor as 
Sij = PYij + Tij (21) 


where the second term is trace-free by construction and represents the anisotropic stress 
component (it is initially non-zero for the Alcubierre metric). We then impose that 

O01; = BE Dio; + apj" + Tip; p" — Kad;;, (22) 

Xp = pap + Ap — B'd;p — kalp — p), (23) 
where the first terms are the advection of the fields in spacetime and the terms containing 
k set a timescale related to proper time for the decay of the anisotropic stresses to zero 
and the pressure to that of a stiff fluid for which p = p. We take the timescale 1/« to 
be of order [R]. Evolving the anisotropic and average pressure components, we can at 
any point reconstruct S;; from their values as defined by Eq. (21)*. 

This construction assumes that there is some transition after the warp field 
breakdown that reverts the warp fluid to a (phantom) perfect fluid with a stiff equation 
of state and a velocity that matches the normal observers, such that the divergence 
term in Eq. (19) and the anisotropic stresses decay away. Further investigations are 
needed to understand the influence of the equation of state of the warp field on the 
gravitational-wave signal, but it appears that the metric behaviour mostly dominates 
the signal. However, the chosen equation of state will strongly affect the matter fluxes, 
which may change the conclusions in Sec. 3. 

The evolution of the energy density p that results from this prescription is illustrated 
in Fig. 3. 


2.5. Numerical setup 


The size of the domain is L = 168 [R] with N = 224 cells on the coarsest grid in the 
direction of travel of the warp ship (the z axis). In the xz and y directions, we use 
the symmetry of the problem to halve the simulation domain so that we only evolve a 


+ One might consider trying to naively evolve the equation of state in Eq. (19), but a simple analysis 
shows that this introduces higher derivative terms to the equation of motion for the momentum density 
(either 3rd spatial derivatives of the lapse or 3rd-time derivatives of the spatial metric). Given the high 
chance that such terms will result in an ill-posed evolution system, we do not follow this approach. 
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quarter of the full domain, with appropriate parities at the reflection boundaries. In the 
nonsymmetric directions, we use Sommerfeld boundary conditions; cf. [5]. We employ 6 
levels of 2:1 refinement, such that the finest grid has a width of 2.625[R] and the finest 
grid spacing is dx ~ 0.0117[R]. 

For all runs, the size of the warp ship in code units is set to be R = 1, and the 
value that determines the wall thickness ø is set to be 8. The velocity of the warp ship 
is varied between 0.1 and 0.2. We have found that we can stably evolve higher velocities 
up to 0.5, but good quality results require very high resolutions, which would be better 
achieved with a code adapted to the cylindrical symmetry. Above v = 1, we expect to 
encounter pathologies related to the ship’s faster-than-light speed, but again, we would 
require higher resolutions and more computational resources to confirm this. 

Details of the checks we have done on convergence and verification of the matter 
evolution are given in Appendix A. 


3. Results 


3.1. Gravitational waves 


A visualisation of the evolution of the gravitational-wave content is shown in Fig. 2. 
We plot the extracted Weyl scalar V4 for the main m = 0 modes in Fig. 4. The 
signal comes as a burst, initially having no gravitational wave content, followed by an 
oscillatory period with a characteristic frequency of order 1/{R]. Overall, the signal is 
very distinct from the typical compact binary coalescences observed by gravitational- 
wave detectors and more similar to events like the collapse of an unstable neutron star 
or the head-on collision of two black holes [58, 59]. However, since the remnant is not 
a black hole, there is no characteristic black hole ringdown in the signal, and there is 
evidence of a longer lower frequency tail, particularly for higher velocities. 

Integrating Y4, we find that for a warp speed of v = 0.1, the magnitude of the 
strain times the radius of extraction is of order 10~*[R]. Thus, for a 1 km-sized warp 
ship, the strain of the signal at a distance of 1 Mpc would be 10~7!. However, its 
frequency (f ~ 1/[R]) equates to f ~ 300 kHz, which is outside the range of current 
ground based detectors. Proposals for higher frequency detectors have been made [60], 
so in future one may be able to put bounds on the existence of such signals. Depending 
on the size of the warp bubble, which determines the frequency of the signal, as well as 
the distance to the source, these signatures may be detected through a coherent search 
or burst search combining information from multiple gravitational-wave detectors, e.g., 
[61-63]. A more in depth analysis, like the one of [64], would be needed to give accurate 
bounds. 


3.2. Matter fluxes 


A visualisation of the evolution of the matter-energy content is shown in Fig. 3. In 
Fig. 5, we show a comparison of the matter and gravitational-wave fluxes out of the 
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t= Vert 


Figure 4. Modes of the Newman-Penrose scalar Y4 multiplied by the coordinate 
radius for different values of the warp speed v. Given the symmetry of the problem, 
the m4 0 modes are zero, so we only show the modes with l = 2, 3, 4. 


sphere at a coordinate radius of r = 14[R] around the warp bubble versus time. We 
see that whilst the gravitational-wave flux is strictly positive as expected, the matter 
flux alternates in sign, with waves of positive and negative energy leaving the volume. 
It appears from the evolution that the ejection of matter is driven by a shell of positive 
energy forming within the negative energy density shell of the initial bubble and driving 
it outwards. This is reminiscent of the early works of Bondi [42], showing the effect of 
repulsion and attraction between positive and negative energy matter. 

In Fig. 6, we show a zoom-in of the fluxes and integrated energy loss, which shows 
that the matter flux results in an overall negative flux out of the spacetime, with the 
gravitational waves giving a much smaller positive flux. The result is that the overall 
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Figure 5. A comparison of the matter and gravitational-wave energy fluxes out of 
sphere at a coordinate radius of r = 14[R] around the warp bubble versus time for 
different v. The left panel shows the integrated GW fluxes in the 20 mode, whereas 
the right shows the integrated matter fluxes. The net gravitational-wave flux is always 
positive, whereas the net matter flux oscillates as waves of positive and negative energy 
leave the volume. 
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Figure 6. Zoom in on the start and end periods of the integrated energy fluxes in the 
matter and gravitational waves as in Fig. 5 for the case v = 0.1. We see that the flux 
begins in the matter sector, with alternating waves of positive and negative energy. 
The gravitational-wave flux starts a little later but continues throughout the periods 
in which the net matter flux is both positive and negative. The second panel shows 
that the overall matter energy loss is negative (i.e. the energy of matter in the volume 
increases), whereas the gravitational-wave flux is overall positive as expected but much 
smaller. The final ADM mass of the volume is therefore larger than its initial value 
(which was zero). 


ADM mass of the spacetime will become more positive (having started at zero). It 
would be interesting to study whether this is a generic effect, such that the ADM 
energy for spacetimes with NEC-violating matter does not dynamically decrease below 
zero. However, we can see from Fig. 6 that nothing appears to impede the gravitational- 
wave flux when the net matter flux is instantaneously positive overall, as between times 
t — Text = —1.0 and —0.5 on the left panel of the plot. 

Since we do not know the type of matter used to construct the warp ship, we do not 
know whether it would interact (apart from gravitationally) with normal matter as it 
propagates through the Universe. If it did interact, it may give rise to further signatures 
(i.e., a multimessenger event). For a 1 km sized warp bubble traveling at v = 0.1, the 
magnitude of the energy carried by the matter waves, of order 10~? [R] (see Fig. 5), 
would be around 1/100 times the mass-energy of the sun. 
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4. Discussion 


In this work, we proposed a formalism for studying warp drive spacetimes dynamically 
and produced the first fully consistent numerical-relativity waveforms for the collapse 
of a warp drive bubble. 

We also studied the evolution of the fluxes of matter and gravitational-wave energy 
from the spacetime, finding that an initial wave of negative energy matter is emitted, 
followed by alternating positive and negative energy waves. The gravitational-wave 
flux starts shortly after and is positive throughout, as expected. The end result is a 
net negative flux of energy, resulting in an increased ADM mass of the spacetime. It 
would be interesting to explore to what extent this behaviour is generic or depends on 
the chosen equation of state supporting the warp bubble. To do this, one would want 
to consistently implement a fluid description with a general rest frame velocity and 
equation of state, unlike the somewhat ad hoc setup employed here. 

This work could be extended in several other directions. Firstly, using a code better 
adapted to the cylindrical symmetry, one could extend to higher speeds, including 
those that exceed v = 1.0, to see if there is a change in behaviour as light speed is 
approached. One might expect to see some kind of pathology develop, e.g., the formation 
of elliptic regions that cause the dynamical evolution to fail. Better resolution and longer 
evolutions would also allow better characterization of the strain signal, particularly the 
role of the gravitational-wave memory content (see [65]). 

As discussed above, for a lkm-sized ship, the frequency of the signal is much 
higher than the range probed by existing detectors, and so current observations cannot 
constrain the occurrence of such events. However, the amplitude of the strain signal 
would be significant for any such event within our galaxy and even beyond, and so 
within the reach of future detectors targeting higher frequencies [60]. We caution 
that the waveforms obtained are likely to be highly specific to the model employed, 
which has several known theoretical problems, as discussed in the Introduction. Further 
work would be required to understand how generic the signatures are, and properly 
characterise their detectability. 
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Appendix A. Numerical details, validation and convergence 


To validate our simulations, we monitor the violations of the Hamiltonian and 
momentum constraints which are given by 


H:=R+ SK? — Ay A™ — 2A — 167p, (Ast) 
ieee 7 2 
X 


In the continuum limit these two quantities should vanish. We check that the initial 
values converge to zero at the appropriate order, which validates the initial stress energy 
tensor components p and S;. The Hamiltonian constraint is solved exactly to numerical 
precision, whereas the spatial gradients in the Momentum constraint mean that it should 
initial converge at 4th order. During the evolution we expect errors to be dominated 
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by the time integration, and give something between 1st and 3rd order, with some of 
the drop in convergence caused by the approximate way that we zero the constraints 
outside the volume. 
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Figure A1. The top panel shows the L2 norm of the Hamiltonian constraint violation 
over time in the central volume over which data is extracted, the second the same 
for the Momentum constraint violation. These show approximate convergence of the 
quantities to zero with increasing resolution between Ist and and 3rd order. The 
final panel shows the reconciliation of the matter-energy fluxes and source term (the 
exchange of energy with the curvature) as described in the main text. The agreement 
between the direct measurement of M and its inferred value from the flux and source 
terms validates that the matter is well resolved in the relevant volume. 


We calculate the L2 norm of the constraints over the coordinate volume for which 
the matter integrals are evaluated (with its 2D surface also being the surface used for 
extraction of the gravitational waves), which is at a coordinate radius of r = 14[R] from 
the centre of the grid. As shown in Fig Al, the constraint violations remain under 
control throughout the simulation and converge at a reasonable order. 

To verify that the matter is sufficiently resolved within the volume of interest, we 
also reconcile the matter quantities over time, checking that the net energy flux out 
of the surface matches the change in the matter energy and any exchange with the 
curvature (often referred to as a “source”) over time. For this, we follow the approach 
of [66]. As shown in the last panel of Fig. Al, the agreement is good throughout the 
simulation, which gives us confidence in the matter evolution. 


References 


[1] Abbott B P et al. (LIGO Scientific, Virgo) 2018 (Preprint 1811.12907) 
[2] Abbott R et al. (LIGO Scientific, Virgo) 2021 Phys. Rev. X 11 021053 (Preprint 
2010.14527) 


[3] Abbott R et al. (KAGRA, VIRGO, LIGO Scientific) 2023 Phys. Rev. X 13 041039 
(Preprint 2111.03606) 


REFERENCES 15 


4] Gourgoulhon E 2007 arXiv: 0703035 (Preprint gr-qc/0703035) 


5] Alcubierre M 2008 Introduction to 3+1 Numerical Relativity (Oxford University 
Press) 


6] Baumgarte T and Shapiro S 2010 Numerical Relativity (Cambridge: Cambridge 
University Press) 


7| Baumgarte T W and Shapiro S L 2021 Numerical Relativity: Starting from Scratch 
(Cambridge University Press) ISBN 978-1-108-93344-5, 978-1-108-84411-6, 978-1- 
108-92825-0 


[8] Chandrasekhar S 1931 Astrophys. J. 74 81-82 
[9] Penrose R 1965 Phys. Rev. Lett. 14 57-59 


10] Akiyama K et al. (Event Horizon Telescope) 2019 Astrophys. J. Lett. 875 L1 
(Preprint 1906. 11238) 


11] Abbott B P et al. (LIGO Scientific, Virgo) 2016 Phys. Rev. Lett. 116 061102 
(Preprint 1602 .03837) 


12] Ghez A M et al. 2008 Astrophys. J. 689 1044-1062 (Preprint 0808 . 2870) 


13] Lobo F S 2007 Exotic solutions in General Relativity: Traversable wormholes and ’warp drive’ spacet 
(Preprint 0710.4474) 


14] Santiago J, Schuster S and Visser M 2021 Universe 7 271 (Preprint 2106 .05002) 
15] Alcubierre M 1994 Class. Quant. Grav. 11 L73-L77 (Preprint gr-qc/0009013) 
16] Natário J 2002 Class. Quant. Grav. 19 1157-1166 (Preprint gr-qc/0110086) 


17] Van Den Broeck C 1999 Class. Quant. Grav. 16 3973-3979 (Preprint gr-qc/ 
9905084) 


18] Krasnikov S 1998 Phys. Rev. D 57 4760-4766 (Preprint gr-qc/9511068) 

19] Everett A E and Roman T A 1997 Phys. Rev. D 56 2100-2108 (Preprint 
gr-qc/9702049) 

20) Alcubierre M and Lobo F S N 2017 Fundam. Theor. Phys. 189 257-279 

21] Everett A E 1996 Phys. Rev. D 53 7365-7368 

22| Hiscock W A 1997 Class. Quant. Grav. 14 L183-L188 (Preprint gr-qc/9707024) 
23| Visser M, Bassett B and Liberati S 2000 Nucl. Phys. B Proc. Suppl. 88 267-270 
(Preprint gr-qc/9810026) 

24| Loup F, Waite D and Halerewicz E J 2001 (Preprint gr-qc/0107097) 

25] Obousy R K and Cleaver G 2008 J. Br. Interplanet. Soc. 61 364 (Preprint 
0712.1649) 

26] Olum K D 1998 Phys. Rev. Lett. 81 3567-3570 (Preprint gr-qc/9805003) 

27| Lobo F S N and Visser M 2004 Class. Quant. Grav. 21 5871-5892 (Preprint 
gr-qc/0406083) 

28] Lentz E W 2021 Class. Quant. Grav. 38 075015 (Preprint 2006 .07125) 


REFERENCES 16 


29 


30 


31 


[32 


33 


34 


35 
36 
37 


38 


39 


40 
41 
42 
[43 
44 


45 


46 


47 


48 


49 


50 


Bobrick A and Martire G 2021 Class. Quant. Grav. 38 105009 (Preprint 2102. 
06824) 

Fell S D B and Heisenberg L 2021 Class. Quant. Grav. 38 155020 (Preprint 
2104. 06488) 

Santiago J, Schuster S and Visser M 2022 Phys. Rev. D 105 064038 (Preprint 
2105 . 03079) 

Barcelo C and Visser M 2002 Int. J. Mod. Phys. D 11 1553-1560 (Preprint 
gr-qc/0205066) 

Kontou E A and Sanders K 2020 Class. Quant. Grav. 37 193001 (Preprint 
2003.01815) 

Pfenning M J and Ford L 1997 Class. Quant. Grav. 14 1743-1751 (Preprint 
gr-qc/9702026) 

Fewster C J and Teo E 1999 Phys. Rev. D 59 104016 (Preprint gr-qc/9812032) 

Krasnikov S 2003 Phys. Rev. D 67 104013 (Preprint gr-qc/0207057) 

Finazzi S, Liberati S and Barcelo C 2009 Phys. Rev. D 79 124017 (Preprint 
0904.0141) 

Alcubierre M 2017 Wormholes, Warp Drives and Energy Conditions vol 189 
(Springer) (Preprint 2103 .05610) 

Schuster S, Santiago J and Visser M 2023 Gen. Rel. Grav. 55 14 (Preprint 
2205 . 15950) 

Witten E 1981 Commun. Math. Phys. 80 381 

Schon R and Yau S T 1979 Commun. Math. Phys. 65 45-76 

Bondi H 1957 Rev. Mod. Phys. 29 423-428 

Arnotwitt R, Deser S and Misner C 1960 Physical Review 117 1595-1602 


York J W 1979 Kinematics and dynamics of general relativity Sources of 


gravitational radiation ed Smarr L L (Cambridge, UK: Cambridge University Press) 
pp 83-126 ISBN 0-521-22778-X 

Drexler D and Clark M 2014 Ships of the Line Star Trek (Gallery Books) ISBN 
9781476792675 URL https: //books. google.co.uk/books?id=9YNOAWAAQBAJ 
Alic D, Bona-Casas C, Bona C, Rezzolla L and Palenzuela C 2012 Phys.Rev. D85 
064040 (Preprint 1106.2254) 

Bona C, Masso J, Seidel E and Stela J 1995 Phys. Rev. Lett. 75 600-603 (Preprint 
gr-qc/9412071) 

Baker J G, Centrella J, Choi D I, Koppitz M and van Meter J 2006 Phys. Rev. 
Lett. 96 111102 (Preprint gr-qc/0511103) 

Campanelli M, Lousto C O, Marronetti P and Zlochower Y 2006 Phys. Rev. Lett. 
96 111101 (Preprint gr-qc/0511048) 

van Meter J R, Baker J G, Koppitz M and Choi D I 2006 Phys. Rev. D73 124011 
(Preprint gr-qc/0605030) 


REFERENCES 17 


[51] 
[52] 


[53] 


54 
59 
56 


57 
58 
59 


60 
61 
62 


63 
64 
65 
66 


Andrade T et al. 2021 J. Open Source Softw. 6 3703 (Preprint 2201 . 03458) 


Radia M, Sperhake U, Drew A, Clough K, Figueras P, Lim E A, Ripley J L, 
Aurrekoetxea J C, França T and Helfer T 2022 Class. Quant. Grav. 39 135006 
(Preprint 2112.10567) 


Clough K, Figueras P, Finkel H, Kunesch M, Lim E A and Tunyasuvunakool S 
2015 Class. Quant. Grav. 32 245011 [Class. Quant. Grav.32,24(2015)] (Preprint 
1503. 03436) 


Alcubierre M 1997 Phys. Rev. D 55 5981-5991 (Preprint gr-qc/9609015) 
Alcubierre M 2003 Class. Quant. Grav. 20 607-624 (Preprint gr-qc/0210050) 


Baumgarte T W and Hilditch D 2022 Phys. Rev. D 106 044014 (Preprint 
2207 . 06376) 


Baker J G, Campanelli M and Lousto C O 2002 Phys. Rev. D 65 044001 (Preprint 
gr-qc/0104063) 


Giacomazzo B, Rezzolla L and Stergioulas N 2011 Phys.Rev. D84 024022 (Preprint 
1105.0122) 


Sperhake U, Kelly B J, Laguna P, Smith K L and Schnetter E 2005 Phys. Rev. D 
71 124042 (Preprint gr-qc/0503071) 


Aggarwal N et al. 2021 Living Rev. Rel. 24 4 (Preprint 2011.12414) 
Sutton P J et al. 2010 New J. Phys. 12 053034 (Preprint 0908.3665) 


Cornish N J and Littenberg T B 2015 Class. Quant. Grav. 32 135012 (Preprint 
1410.3835) 


Klimenko S et al. 2016 Phys. Rev. D 93 042004 (Preprint 1511.05999) 
Kuwahara S and Cannon K 2023 (Preprint 2310. 16067) 

Mitman K et al. 2024 (Preprint 2405 .08868) 

Clough K 2021 Class. Quant. Grav. 38 167001 (Preprint 2104.13420) 


